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Water-leaving  radiances,  retrieved  from  in  situ  or  satellite  measurements,  need  to  be  corrected  for  the 
bidirectional  properties  of  the  measured  light  in  order  to  standardize  the  data  and  make  them  compar¬ 
able  with  each  other.  The  current  operational  algorithm  for  the  correction  of  bidirectional  effects  fiem  the 
satellite  ocean  color  data  is  optimized  for  typical  oceanic  waters.  However,  versions  of  bidirectional  re¬ 
flectance  correction  algorithms  specifically  tuned  for  t3^ical  coeistal  waters  and  other  case  2  conditions 
are  particularly  needed  to  improve  the  overall  quality  of  those  data.  In  order  to  anal3rze  the  bidirectional 
reflectance  distribution  function  (BRDF)  of  case  2  waters,  a  dataset  of  typical  remote  sensing  reflectances 
was  generated  through  radiative  transfer  simulations  for  a  large  range  of  viewing  and  illumination  geo¬ 
metries.  Based  on  this  simulated  dataset,  a  case  2  water  focused  remote  sensing  reflectance  model  is 
proposed  to  correct  above-water  and  satellite  water-leaving  radiance  data  for  bidirectional  effects. 
The  proposed  model  is  first  validated  with  a  one  year  time  series  of  in  situ  above-water  measurements 
acquired  by  collocated  multispectral  and  hyperspectral  radiometers,  which  have  different  viewing  geo¬ 
metries  installed  at  the  Long  Island  Sound  Coastal  Observatory  (LISCO).  Match-ups  and  intefcompar- 
isons  performed  on  these  concurrent  measurements  show  that  the  proposed  algorithm  outperforms  the 
algorithm  currently  in  use  at  all  wavelengths,  with  average  improvement  of  2.4%  over  the  spectral  range. 
LISCO’s  time  series  data  have  also  been  used  to  evaluate  improvements  in  match-up  comparisons  of 
Moderate  Resolution  Imaging  Spectroradiometer  satellite  data  when  the  proposed  BRDF  correction 
is  used  in  lieu  of  the  current  algorithm.  It  is  shown  that  the  discrepancies  between  coincident  in-situ 
sea-based  and  satellite  data  decreased  by  3.15%  with  the  use  of  the  proposed  algorithm.  This  confirms 
the  advantages  of  the  proposed  model  over  the  current  one,  demonstrating  the  need  for  a  specific  case  2 
water  BRDF  correction  algorithm  as  well  as  the  feasibility  of  enhancing  performance  of  current  and  fu¬ 
ture  satellite  ocean  color  remote  sensing  missions  for  monitoring  of  typical  coastal  waters. 


1.  Introduction 

The  water-leaving  radiance  field  emerging  from  sea 
water  is  not  generally  isotropic  because  of  the  noni¬ 
sotropic  character  of  the  volume  scattering  function 
(VSF)  of  the  water  body  coupled  with  nonisotropic 
illumination  conditions  that  generally  prevail  above 
the  sea  surface.  Furthermore,  the  wavy  nature  of  the 
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water  surface  also  has  noticeable  influence  over  the 
anisotropy  of  the  water-leaving  radiance  field.  The 
bidirectional  structure  of  the  emerging  radiance  field 
has  implications,  particularly  in  the  processing  of 
data  from  ocean  color  satellite  sensors  as  well  as 
in  calibrationA^alidation  activities  for  ocean  color 
sensors  in  which  water-leaving  radiances  obtained 
under  various  geometrical  configurations  have  to 
be  quantitatively  and  accurately  compared. 

Variations  in  the  bidirectional  properties  of  the  up¬ 
ward  radiative  field  have  already  been  extensively 
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studied.  Morel  and  Gentili  U]  reported  the  variations 
in  the  bidirectional  properties  of  the  upward  radia¬ 
tive  field  for  case  1  waters,  and  their  studies  related 
the  bidirectional  variations  of  the  water-leaving  ra¬ 
diance  field  to  chlorophyll  concentration,  [Chi],  Their 
findings  were  also  corroborated  by  field  measure¬ 
ments  obtained  using  an  underwater  radiance  distri¬ 
bution  camera  system  in  low-chlorophyll  case  1 
waters  [2].  The  anisotropy  of  the  water-leaving  radi¬ 
ance  was  also  discerned  from  the  series  of  BRDF 
measurements  carried  out  from  research  aircraft 
for  both  case  1  and  case  2  waters  [3].  These  airborne 
measurement  results  show  good  overall  agreement 
with  theoretical  simulations.  Also  based  on  the  field 
measurements  performed  exclusively  for  case  1 
water  types,  Voss  and  Morel  reported  that  the  Q  fac¬ 
tor,  defined  as  the  ratio  between  upwelling  irradi- 
ance  and  upwelling  radiance  from  the  sea  surface, 
can  vary  within  7%  (41.  These  predictive  investiga¬ 
tions  were  tentatively  extended  to  typical  case  2 
waters,  usually  dominated  by  sediment  or  by  colored 
dissolved  organic  matter  (CDOM),  and  are  currently 
being  used  in  the  NASA  ocean  color  data  processing 
procedures  as  the  operational  correction  algorithm 

[5]  for  the  retrieval  of  water-leaving  radiance  from 
both  case  1  and  case  2  water  types.  Although  the  cur¬ 
rent  operational  algorithm  is  considered  reasonably 
effective  for  case  1  water  conditions,  the  validity  of 
this  correction  algorithm  in  case  2  water  types  has 
not  been  really  verified  in  natural  environments 

[6] .  Accordingly,  the  need  for  an  improved  version 
of  a  bidirectional  reflectance  correction  algorithm 
particularly  tuned  for  typical  case  2  water  conditions 
is  the  general  consensus  among  the  ocean  color  radio- 
metry  community  [6,7]. 

In  case  2  waters,  the  more  complex  variability  of 
the  water  composition  results  in  the  bidirectional 
properties  of  the  upwelling  radiance  field  being  se¬ 
verely  affected  by  concentrations  of  CDOM  and 
mineral  particles,  which  can  have  significant  conse¬ 
quences  on  the  correction  factors  both  spectrally  and 
angularly  [71.  Consequently,  extensive  studies  of  the 
bidirectional  properties  of  case  2  waters  have  been 
made  through  radiative  transfer  simulations  with  in¬ 
puts  of  inherent  optical  properties  (lOPs)  typical  for 
case  2  water  types.  Lee  et  al.  proposed  a  semianaly- 
tical  model  that  relates  nadir  viewing  remote  sensing 
reflectance  to  single  backscattering  albedo  and  the 
solar  zenith  angle  for  both  deep  and  shallow  waters 
[8,9).  This  model  was  later  extended  as  both  viewing 
and  solar  zenith  angle  dependent  to  be  used  for  mea¬ 
surements  made  at  various  viewing  angles  [9).  Para¬ 
meterization  of  the  remote  sensing  reflectance  as  a 
function  of  single  backscattering  albedo,  solar  zenith 
and  sensor  viewing  angles,  and  the  sea  surface  wind 
speed  was  also  developed  by  Albert  and  Mobley  [10] 
and  Albert  and  Gege  [111.  Nevertheless,  these  mod¬ 
els  do  not  take  into  account  the  solar-sensor  relative 
azimuth  angle  as  an  input  parameter,  and  as  a  re¬ 
sult,  bidirectional  variations  in  the  azimuthal  direc¬ 
tion  cannot  be  analyzed  with  the  use  of  these  models. 


Notably,  Park  and  Ruddick  [6]  proposed  a  remote 
sensing  reflectance  model  based  on  lookup  tables 
of  coefficients  generated  for  a  large  range  of  solar  ze¬ 
nith,  sensor  viewing,  and  solar-sensor  relative  azi¬ 
muth  angles  to  estimate  the  directional  variation 
of  the  water-leaving  radiance  field  of  both  case  1 
and  case  2  waters.  However,  their  model  requires  a 
phase  function  parameter  defined  by  the  contribu¬ 
tion  of  suspended  particles  to  the  backscattering 
coefficients.  But  in  remote  sensing  of  coastal  waters, 
accurate  retrieval  of  such  parameters  of  the  water 
body  remains  highly  challenging,  and  this  can  limit 
the  operational  use  of  the  algorithm.  Lee  et  al.  [12] 
also  proposed  an  approach  centered  on  the  lOP  of 
water  to  remove  the  angular  Variations  in  remote 
sensing  reflectance  measurements.  The  evaluations 
of  the  overall  performance  of  their  BRDF  correction 
scheme  carried  out  with  the  use  of  the  field  measured 
data  indicate  slight  improvement  over  the  current 
standard  correction  algorithm. 

In  our  previous  studies  [7],  bidirectional  reflec¬ 
tance  effects  of  coastal  waters  were  anal3^ed  through 
extensive  radiative  transfer  simulations  of  remote 
sensing  reflectances  t5q)ical  for  a  wide  range  of  coast¬ 
al  water  conditions.  Based  oti  this  analysis,  we 
proposed  a  model  that  relates  the  single  backscatter¬ 
ing  albedo  to  the  remote  sensing  reflectance  through 
a  leaist  mean  squared  optimization  for  a  large  range 
of  lOPs  and  of  viewing  and  illumination  geometries. 
Furthermore,  the  bidirectional  reflectance  param¬ 
eter  Q  was  also  analyzed  based  on  a  microphysical 
model  of  the  water  constituents  and  was  observed 
to  have  variations  with  a  wavelength  of  about  15%, 
esj>ecially  in  the  blue  and  green  regions  of  the  water¬ 
leaving  radiance  spectrum,  depending  on  the  water 
composition. 

As  an  extension  of  the  previous  model,  we  propose, 
in  this  paper,  an  improved  version  of  the  remote  sen¬ 
sing  reflectance  model  based  on  lookup  tables  gener¬ 
ated  for  a  large  set  of  illumination  and  viewing 
geometries  and  water  optical  properties,  to  estimate 
the  directional  variation  of  the  water-leaving  radi¬ 
ance  field  of  t)q)ical  case  2  water  conditions.  Unlike 
previous  studies  [6,7,12],  the  model  coefficients  were 
generated  as  being  wavelength  dependent  in  order  to 
better  deal  with  the  spectral  dependency  of  bidirec¬ 
tional  effects.  The  proposed  model  is  applied  to  two 
different  above-water  datasets  acquired  by  collo¬ 
cated  instruments  on  the  Long  Island  Sound  Coastal 
Observatory  (LISCO).  The  LISCO  site,  near  North- 
port,  New  York,  has  been  established  to  support  pre¬ 
sent  and  future  multispectral  and  h)q)erspectral 
ocean  color  sensor  calibrationAralidation  activities, 
as  well  as  the  development  of  new  measurement 
and  retrieval  techniques  and  algorithms  for  coastal 
waters  [13-15].  These  sea-based  datasets  also  permit 
us  to  examine  and  validate  satdlite  retrievals  based 
on  the  proposed  and  the  current  BRDF  algorithms. 

In  Section  2,  fundamental  equations  and  rationale 
of  remote  sensing  reflectance  retrieval  are  briefly  re¬ 
visited.  The  following  section,  Section  3,  deals  with 
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the  radiative  transfer  modeling  of  coastal  waters  and 
the  details  of  the  proposed  BRDF  correction  scheme. 
Theoretical  and  noise  sensitivity  analysis  based  on 
the  simulated  dataset  is  described  in  Section  4. 
The  performances  of  the  new  scheme  are  tested  in 
Section  5  based  on  time  series  of  coincident  remote 
sensing  reflectance  measurements  at  different  scat¬ 
tering  angles  from  the  LISCO  platform.  Finally,  the 
proposed  BRDF  correction  scheme  is  applied  to 
Moderate  Resolution  Imaging  Spectroradiometer 
(MODIS)  satellite  data,  and  its  performance  is  dis¬ 
cussed  based  on  comparisons  with  the  sea-based 
LISCO  datasets.  This  is  followed  by  concluding 
remarks. 


2.  Theoretical  Background 

The  fundamental  equation,  shown  below,  relates  any 
slanted  radiance  exiting  the  water,  L^,,  to  the  down- 
welling  irradiance  just  above  the  sea  surface,  and 
to  the  absorption  and  backscattering  coefficients  of 
the  water  body  (a  and  6^)  through  three  quantities: 
31  (a  dimensionless  factor  that  merges  all  the  reflec¬ 
tion  and  refraction  effects  that  occur  when  downward 
irradiance  and  upward  radiance  propagate  through 
the  air-water  interface),  f  (a  dimensionless  coeffi¬ 
cient  that  relates  the  magnitude  of  the  irradiance  re¬ 
flectance  just  below  the  surface  to  the  lOPs,  namely  a 
and  bfy)  and  Q  (the  bidirectional  function  with  vmits  of 
sr"^)  as  follows  [5]: 


Lu,(Os. A,  w,  lOP)  =  X)  x  m{0,, 

f{0,,X,W^OP) 

X  .  (1) 

a  (A)  -h 

The  water-leaving  radiance,  at  a  given  wave¬ 
length,  Xy  depends  on  its  own  direction,  depicted  by 
Of)y  the  viewing  zenith  angle,  and  <py  the  relative  azi¬ 
muth  angle  from  the  Sim's  direction,  on  the  illumina¬ 
tion  conditions  (determined  by  the  solar  zenith  angle, 
Og)y  on  the  sea  state  (parameterized  via  the  wind 
speed,  W)y  and  on  the  lOPs  of  the  water  body  itself 
In  contrast,  besides  atmospheric  parameters,  is 
only  a  function  of  the  solar  zenith  angle,  Ogy  and 
the  wavelength.  The  last  factor,  6^,/(a  -h  6^),  is  custo¬ 
marily  called  single  backscattering  albedo,  and  it  is 
routinely  denoted  as  co.  In  Eq.  (^)  is  the  in-water 
viewing  angle,  and  it  is  related  to  0^  through  Snell’s 
law  (sin  6^  =  n  sin  where  n  is  the  refractive  index 
of  water). 

The  above-water  remote  sensing  reflectance  (Rrs) 
is  defined  as  [161 


RT3{0gy0^y<f)yXyWjOP) 


lOP) 

Ed(PayX) 


(2) 


Therefore,  from  Eqs.  (1.)  and  (2),  Rrs  can  be  related  to 
fy  Qy  and  the  lOPs  of  the  water  as  follows: 


mgyXyWylOP) 

""  QiOgMy^^y^ylOV) 

X  (I}{X).  (3) 

On  the  other  hand,  Gordon  et  al.  [17]  showed  that 
nadir  viewed  remote  sensing  reflectance  can  be  mod¬ 
eled  as  the  polynomial  flmctioti  of  cd.  It  was  also 
shown  in  [10,1^]  that  the  factor//Q  can  be  parame¬ 
terized  as  the  function  of  tw,  Ogy  and  for  both  deep 
and  shallow  waters.  In  fact,  cd  Can  be  considered  as 
the  representative  of  the  lOP  input  argument  to 
the  f/Q  factor,  and  the  relationship  between  them 
has  to  be  subsequently  expressed  in  linear  or  higher 
order  polynomial  functions  [7,jJ0,j^].  Consequently, 
the  relationship  between  Rrs(A)  and  a){X)  is  also  con¬ 
sidered  to  be  a  second  or  higher  degree  j>ol3momial 
[6,7,10,12,17-19],  In  this  study,  we  model  Rrs(A)  as 
the  third  order  pol3nnomial  of  g>(A)  as  follows: 

3 

1  =  1 


and  the  model  coefficients,  a,,  of  the  relationship 
[Eq.  (4)1  are  tabulated  for  a  large  range  of  solar  ze¬ 
nith  (0g)y  viewing (^j^),  and  relative  azimuth  i<p)  angles 
as  well  as  for  wavelengths  (A).  This  poljmomial  model 
merges  all  bidirectional  effects  and  directly  relates 
the  Rrs  to  the  single  parameter  allowing  the  BRDF 
correction  to  be  performed  using  only  the  informa¬ 
tion  on  the  viewing  and  illumination  geometries  at 
the  time  of  the  Rrs  measurements.  Thus,  the  model 
offers  a  computationally  straightforward  and  simple 
approach  for  performing  the  directional  correction  by 
forgoing  the  potential  to  model  the  different  influ¬ 
ences  onLy;  (and  consequently  on  the  Rrs)  separately 
as  in  Eqs.  (D  and  (3). 

Based  on  the  simulated  dataset,  which  will  be  de¬ 
tailed  in  the  following  section,  it  was  found  that  high¬ 
er  order  polynomials  render  les$  least  square  (LSQ) 
curve  fitting  errors.  Nevertheless,  it  was  also  ob¬ 
served  that  higher  than  third  order  pol3momial  fit¬ 
tings  also  show  tendencies  to  give  unrealistic  Rrs 
values  if  the  input  co  values  were  outside  the  range 
of  the  simulated  dataset.  Thus,  we  decided  to  use  the 
third  order  as  the  optimal  pol5nnomial  degree  in  our 
modeling.  The  fundamental  difference  between  the 
proposed  model  and  the  previous  models  is  that,  un¬ 
like  previous  studies,  model  coefficients  of  the  Rrs-ty 
relationship  are  generated  as  wavelength  dependent 
in  order  to  better  deal  with  the  known  spectral  var¬ 
iation  in  the  Q  factor  [7,18].  Thus,  this  proposed  mod¬ 
el  allows  the  BRDF  correction  to  be  performed  for 
each  wavelength  independently,  demonstrating  the 
potential  for  better  dealing  with  the  spectral  depen¬ 
dency  of  the  BRDF  [5,7,18],  In  the  proposed  model, 
the  factor  31  is  integrated  in  the  formulation  by  di¬ 
rectly  relating  o)  to  the  above-water  simulated  Rrs. 
Morel  et  al.  previously  reported  large  variation  in 
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with  wind  speed  [5,M].  Nevertheless,  Gk)rdon^s  re¬ 
vision  [^]  of  their  results  later  showed  that  the  in¬ 
fluence  of  surface  roughness  on  the  transfer  of 
subsurface  radiance  through  the  air-water  interface 
is  minimal  for  surface  viewing  angles  ^  60°.  Albert 
and  Mobley  have  also  studied  the  dependence  of  ir- 
radiance  reflectance  for  infinitely  deep  waters  on  sur¬ 
face  wind,  obtaining  the  same  results  [10],  Also,  Park 
and  Ruddick  [6],  based  on  the  Hydrolight  simula¬ 
tions  for  lOPs  of  both  case  1  and  2  waters  with  the 
different  input  wind  speeds  (0,  5  and  10  ms"^), 
showed  the  effect  of  wind  speed  on  bidirectional  re¬ 
flectance  to  be  small.  Following  these  conclusions,  we 
did  not  incorporate  the  wind  speed  (W)  as  an  input 
parameter,  because  its  influence  on  both  the  transfer 
of  subsurface  radiance  through  the  air-water  inter¬ 
face  and  the  underwater  radiance  field  is  minimal 
112,^].  It  should  be  noted  that  the  BRDF  correction 
usually  accounts  for  10%— 20%  of  Rrs(A)  16,7].  Thus, 
the  estimation  of  the  correction  quality  is  a  challen¬ 
ging  task,  especially  for  the  field  data,  where 
other  factors  can  contribute  significantly  to  total 
uncertainties. 

3.  Coastal  Waters  Radiative  Transfer  Modeling 
3.1.  iOP  Model 

A  dataset  of  remote  sensing  reflectance  spectra, 
Rrs(A),  typical  for  coastal  water  conditions  has  been 
generated  by  us  through  radiative  transfer  simula¬ 
tions  using  Hydrolight  5  [22].  In  these  simulations, 
a  four-component  bio-optical  model  was  assumed. 
The  four  components  are  pure  sea  water,  CDOM, 
nonalgal  particulates  (NAP),  and  chlorophyll - 
containing  particles  (phytoplankton).  We  use  the 
subscripts  “u;,”  “y,”  “NAP,”  and  “Chi”  to  identify  these 
four  components  in  the  text  below.  The  input  para¬ 
meters  for  the  modeled  lOPs  of  the  components 
(except  for  pure  sea  water)  are  generated  as 
uniformly  distributed  random  variables  in  the  pre¬ 
scribed  ranges  typical  for  coastal  water  conditions. 
The  most  important  relationships  used  to  compute 
the  input  IOP  parameters  are  detailed  in  the  follow¬ 
ing  paragraphs. 

The  total  spectral  absorption  coefficient  a(A)  is  gi¬ 
ven  as 

a  (A)  =  a^{k)  -f  a^(A)  +  achi(^)  +  onap(-^)»  (5) 

where  the  water  absorption  spectrum  a„,(A)  was 
obtained  from  Pope  and  Fry  [23].  The  chlorophyll  ab¬ 
sorption  coefficient  for  case  2  waters  was  considered 
proportional  to  [Chi]  as  is  often  assumed  [^]  and  is 
given  by 

achi(>l)  =  [Chl]-a*jj(A),  (6) 

where  specific  chlorophyll  absorption 

coefficient,  in  m^/mg.  With  a  wide  variety  of  specific 
absorption  spectra  available,  in  our  simulations  the 
set  of  specific  chlorophyll  absorptions  was  taken  from 


Ciotti  et  al.  [25]  as  a  sum  of  specific  absorptions  of 
microplankton  (a^-  )  and  pico^lankton  (a^^^)  with 

different  weighting  mctors,  Sfy  hs  follows: 

«chi w  =  Sf  ■  +  (1  *  Sf)  •  a (7) 

where  =  [0.1, 0.2, 0.3].  However,  a  recent  study 
[26]  showed  that  the  use  of  th$se  values  generally 
leads  to  an  underestimation  of  in  the  red  and 
near-IR  parts  of  the  spectrum.  L/Onsequently,  we 
modified  the  spectra  shape  in  the  650-700  nm  part 
of  the  spectrum  by  setting  values  at  675  nm  to 
0.0142, 0.0156,  and  0.02  m^/mg  for  Sf  =  0.1, 0.2,  and 
0.3,  respectively. 

The  absorption  spectra  of  NAP  and  CDOM  were 
both  modeled  as  having  exponentially  decreasing 
magnitudes  with  respect  to  the  wavelength  and  de¬ 
termined  from  their  base  reference  values  at  400  nm 
as  in  [^,M] .  A  minerogenic  particle  type  is  assumed 
for  NAP,  and  its  absorption  spectra  is  modeled  as 
follows: 

GtNAP(^)  =  C^NAP  ’  '  ®^P(“‘5nap(^  “  400)), 

(8) 

where  aj;j^(400)  is  the  specific  absorption  of  NAP  at 
400  nm  in  m^/g,  and  C^ap  is  their  concentration  in 
g/m^,  CDOM  absorption  was  also  modeled  as  shown 
below: 

ay(k)  =  a^(400)  •  exp(-$^(A  -  400)),  (9) 

where  a_y(400)  is  the  absorption  of  CDOM  at  400  nm. 
In  Eqs.  (8)  and  (9),  Snap  Sy  denote  the  spectral 
slopes  for  the  absorption  spectra  of  NAP  and  (I)DOM, 
respectively. 

In  a  similar  way,  the  total  scattering  coefficient 
was  simulated  as  a  sum  of  three  components: 

6(A)  =  6iy(A)  -j-  6chi(^)  +  ^NAp('^)t  (10) 

where  6(A)  is  the  total  scattering  coefficient  in  m”^. 
The  scattering  of  NAP  is  modeled  as  a  power  law 
function  [29]  as  follows: 

/550\^nap 

6napW  =  Cnap  •  6^(550)  •  (— j  ,  (11) 

where  6jj^(550)  is  the  specific  scattering  of  NAP  at 
550  nm,  in  m^/g. 

The  scattering  coefficient  of  phytoplankton  is  cal¬ 
culated  as  the  difference  between  the  attenuation 
and  the  absorption  coefficients:’ 

^ChlC'^)  =  CchlC'i)  “  ttchiC'l)*  (12) 

where  6chi(^)  and  cchi(^)  are  the  scattering  and 
attenuation  spectra  of  phjftoplankton.  Then  the  at¬ 
tenuation  spectrum  is  modeled  as  a  power  law  func¬ 
tion  [30]  as  follows: 
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ccmW  =  P  ■  [Chlf  62  .  ^5^ j .  (13) 

The  backscattering  coefficient  can  be  broken  down 
as  the  sum  of  the  contributions  of  the  main  scattering 
constituents: 


^6^  =  bhwW  +  ft6Chl  *  ^ChlU)  +  66NAP  *  ^NApWj 

(14) 

where  is  obtained  in  accordance  with  [^],  and 

bbch\  65NAP  ^re  the  backscattering  ratios  for 
chlorophyll  and  NAP,  assumed  to  be  independent 
of  wavelength,  is  calculated  according  to  the 
following  empirical  parameterization  proposed  by 
Morel  et  al.  [5]: 

bbch\  =  0.002  +  {0.01[0.5  -  0.25  logio(Chl)]}.  (15) 

This  allows  us  to  use  a  particle  phase  function  that 
is  able  to  change  progressively  with  [Chi]  value,  in 
such  a  way  that  the  backscattering  efficiency  can  de¬ 
crease  with  increasing  [Chi]  [5].  is  taken  as  a 
fixed  value  of  0.0183  assuming  a  strong  mineral  com¬ 
ponent  in  NAP. 

Similar  to  the  procedures  followed  for  total  absorp¬ 
tion  and  scattering  coefficients,  the  VSF,  denoted  as 
fi{y/)y  is  broken  down  as  the  sum  of  the  contributions 
of  pure  water  NAP  and  chlorophyll- 

containing  particles  ifich\)  can  be  written  into 
the  summation  of  all  three  components: 

fiiw)  =  fiwiw)  +  /^napW  +  ficuiw)-  (16) 

Equation  (16)  can  be  also  written  in  terms  of 
scattering  coefficients  and  scattering  phase  functions 
(SPFs)  of  the  contributing  components  as  follows 

mh 

bp{w)  =  bu}~Pu}{w)  4-  ^NAP^NAP(V')  +  ^Chl^Chl(V')>  (17) 

where  P{y/)  is  the  SPF  of  sea  water.  The  SPF  of  pure 
sea  water,  p^j^,  has  been  well  studied  and,  for  our  si¬ 
mulations,  is  taken  from  [31] .  In  contrast  to  the  SPF 
of  pure  sea  water,  SPFs  of  nonalgal  and  chlorophyllic 
particles  are  not  well  established  and  show  high 
variability  with  the  water  types  because  of  the  diver¬ 
sity  of  size  distributions,  refractive  indices,  and 
shape  anisotropies  of  particles  [6].  The  most  com¬ 
monly  used  SPF  is  Petzold’s  average-particle  phase 
function,  derived  from  eight  VSF  measurements  car¬ 
ried  out  by  Petzold  in  San  Diego  Harbor  [33,34]. 
However,  Petzold's  phase  function  results  in  a  con¬ 
stant  backscattering  ratio  of  0.()183,  thereby  prohi¬ 
biting  application  of  variable  fe^chi  values  [which 
decrease  with  increasing  [Chi];  see  Eq.  (15)]. 

To  meet  this  requirement  for  the  backscattering 
ratio,  we  employed  the  phase  functions  given  by 
the  Fournier  and  Forand  (FF)  anal3rtic  expression 
[35],  which  is  physically  based  on  l^e  theory  and 


has  reasonable  assumptions  about  the  index  of  re¬ 
fraction  and  particle  size  distribution  of  oceanic 
particles  [^].  With  the  use  of  the  FF  analytical  ex¬ 
pression,  phase  functions  with  the  desired  backscat¬ 
tering  ratios  can  be  generated  by  setting  the  real 
index  of  refraction  of  the  particles,  and  the  Junge 
slope  parameter  of  their  assumed  distribution,  fip. 
Thus,  ^chi  used  in  our  simulations  are  generated 
as  FF  phase  functions  by  setting  npcw  to  a  fix  value 
of  1.06  while  varying  the  /ipchi  values  to  achieve  the 
desired  fe^chi  values,  as  described  in  [29,34].  It  should 
be  noted  here  that  the  shape  of  the  obtained  wdth 
the  FF  phase  functions  is  very  similar  to  that  of  Mor¬ 
el  et  al.  [5]  obtained  through  the  T-matrix  method 
(assuming  randomly  oriented  spheroidal  prides 
for  all  [Chi]  levels  used  in  our  simulations),  ^nap  is 
also  computed  through  the  FF  Expression  by  setting 
'^pNAP  und  //pNAP  to  1.1  and  3.58,  respectively,  to  give  a 
backscattering  ratio  of  0.0183  as  in  [34].  Figure  1 
shows  FF  scattering  phase  functions  with  a  backscat¬ 
tering  ratio  of  0.0045  (red)  used  for  algal  particles 
with  [Chi]  =  10  mg/m^,  with  a  backscattering  ratio 
of  0.007  (green)  used  for  algal  particles  with 
[Chi]  =  1  mg/m^,  and  with  a  backscattering  ratio 
of  0.0183  (blue)  used  for  NAP,  together  with  the 
well-known  Petzold’s  average-particle  phase  func¬ 
tion  for  comparison. 

3.2.  Input  Parameters  for  Simulation 
lOP  input  parameters  for  the  simulation  are  gener¬ 
ated  as  independent  random  variables  with  a  uni¬ 
form  distribution  in  the  prescribed  ranges.  These 
selected  ranges  of  the  input  parameters  comprehen¬ 
sively  cover  the  realistic  synoptic  scale  water  condi¬ 
tions  for  typical  coastal  waters  [36].  In  addition,  as 
input  parameters  are  each  generated  independently 
from  one  another,  the  resulting  input  dataset  truly 
represents  the  conditions  of  the  coastal  waters,  in 
which  constituents  of  the  water  body  vary  indepen¬ 
dently.  Table  1  describes  the  ranges  of  the  input 


Particle  Scattering  Phase  Function  P('P) 


Fig.  1.  (Color  online)  SPFs  of  particles  used  in  simulations:  FF 
6^  =  0.0045  (red)  is  the  SPF  with  a  backscattering  ratio  of 
0.0045,  used  for  algal  particles  with  [Chi]  =10  mg/m^,  FF  6*  = 
0.007  (green)  is  used  for  algal  particles  with  [Chi]  =  1  mg/m^, 
FF  6ft  =  0.0183  (blue)  is  used  for  NAP,  and  Petzold  6ft  =  0.0183 
(black  dotted  line)  is  Petzold's  average-particle  phase  function. 
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parameters  used  in  the  simulation.  With  these  as¬ 
sumed  parameters,  total  absorption  of  the  water  at 
412  nm,  a(412),  is  in  the  range  0. 1-2.3  m"^,  the  total 
particulate  backscattering  coefHcient  varies  from 
0.001  to  0.104  m‘^  and  a){X)  values  are  in  the  range 
of  0.005  to  0.3. 

Remote  sensing  reflectance,  Rrs,  water-leaving  ra¬ 
diance,  Ly;,  and  sea  surface  reflectance,  p,  were  com¬ 
puted  with  Hydrolight  with  the  following  criteria 
using  the  lOP  parameters  specified  above: 

6  wavelengths:  412,  443,  491,  551,  668,  and  753  nm 
9  solar  zenith  angles,  0  to  80®  at  10®  intervals 
9  sensor  zenith  angles,  0  to  80®  at  10®  intervals 
13  relative  azimuth  angles,  0  to  180®  at  15® 
intervals 

Wind  speed,  W:  5  m/s  with  cloud  free  skies  for  the 
dataset  of  remote  sensing  reflectance  and  radiance. 

The  three  angles  used  as  simulation  inputs  are  de¬ 
fined  as  follows:  solar  zenith  angle,  0g^  and  sensor  ze¬ 
nith  angle,  Ot,,  are  the  angles  between  the  zenith  line 
pointing  upward  and  the  directions  of  the  Sun  and 
the  sensor,  respectively;  the  relative  azimuth  angle, 
(/>y  is  defined  as  the  angular  difference  between  the 
sensor  azimuth  angle  and  the  solar  azimuth  angle 
(i.e.,  the  relative  azimuth  is  =  0®  when  the  sensor 
is  in  opposition  to  the  Sun  and  (p  =  180®  when  the 
Sun  is  exactly  behind  the  sensorX 

In  our  simulations  we  modeled  the  water  as  a 
homogeneous  and  infinitely  deep  medium  for  the 
sake  of  simplicity;  therefore,  this  modeling  is  applic¬ 
able  only  to  the  optically  deep  waters  where  bottom 
reflectance  is  negligible.  Only  elastic  scattering  was 
considered,  and  Raman  inelastic  scattering  and 
fluorescence  emission  were  not  included.  However, 
exclusion  of  Raman  inelastic  scattering  has  only  a 
negligible  impact  on  the  accuracy  of  the  modeling 
because  of  the  relatively  high  particulate  concentra¬ 
tions  considered  [5].  Chlorophyll  fluorescence  should 
not  have  an  impact  on  the  output  of  water-leaving 
radiance,  because  the  selected  wavelengths  are  far 
from  the  fluorescence  emission  peak  (located  around 
685  nm)  [^1.  The  Hydrolight  standard  quad  layout, 
which  essentially  provides  10®  resolution  in  the  ze- 


Table  1.  Pangea  of  the  lop  Model  Input  Parametera 
Used  In  the  SImulationa 


Parameters 

Range 

IChll 

1  ~  10  mg/m^ 

ay{400) 

0  ~  2  m-‘ 

^NAP 

0.01  ~  2.5  mg/m^ 

0.01  ~  0.02  nm-' 

•Snap 

0.007  ~  0.015  run-' 

aNAp(400) 

0.02-0.1  mVg 

0.5  - 1  mVg 

rchi 

0.1  - 1.6 

rNAP 

0.5-2 

p 

0.1 -0.5 

All  parameters  are  randomly  generated  in  the  prescribed 
ranges  with  a  uniform  distribution. 


nith  (0)  and  15®  in  the  azimuth  ((p)  direction,  was 
used.  Hydrolight’s  built-in  clear  sky  radiance  model 
of  Harrison  and  Coombes  [38]  and  irradiance  model 
of  Gregg  and  Carder  [391  were  employed.  Other  at¬ 
mospheric  conditions  such  as  s0a-level  pressure,  re¬ 
lative  humidity,  horizontal  visibility,  and  ozone 
concentration  are  set  at  default  Values,  which  are  de¬ 
scribed  in  the  Hydrolight  5  technical  documenta¬ 
tion  [22,241. 

3.3.  Simulation  Results  and  Analysts 
In  general,  Rrs(A)  is  strongly  correlated  with  (o{k),  as  is 
known  from  the  previous  studies  [6,10,17].  However, 
the  value  of  Rrs(A)  for  a  given  vqlue  ^£w(A)  does  vary 
with  the  illumination  and  viewing  geometries  (i.e., 
solar  and  sensor  angles)  as  discussed  in  Section  2. 
Figure  2  shows  the  simulated  Rr3(A)  data  as  a  function 
of  (o{X)  for  the  three  conditions  of  Solar  and  viewing  an¬ 
gles  at  412  and  55 1  nm.  It  is  observed  that  the  relation¬ 
ships  between  Rrs  (A)  and  q)(X)  have  a  strong 
dependency  on  solar  angles  at  low  relative  azimuth 
((p)  angles  and  that  the  dependency  becomes  gradually 
weaker  with  increasing  (p  values. 

Using  the  simulated  Rrs(A)  aitd  ft>(A)  data,  we  com¬ 
puted  the  sets  of  tabulated  coefficients  for  the 
gridded  values  of  0^,^  tp^  and  A  according  to  Eq.  (4). 
For  each  grid,  we  used  500  data  points  of  Rrs  and  w  in 
a  third  order  polynomial  least^mean-square  fitting 
procedure  to  produce  the  coefficients  corresponding 
to  the  respective  geometries  and  wavelengths. 
Table  2  shows  the  coefficients  for  the  nadir  viewing 
and  solar  angles  that  are  routinely  employed  in  the 
estimation  of  Rrs(A)  from  (o{X). 

Figure  3  and  Fig.  4  illustrate  the  veiriability  of  the 
coefficients  a,  for  the  whole  range  of  viewing  and  il¬ 
lumination  geometries.  Figure  8  highlights  the  de¬ 
pendency  of  the  coefficients  a,  On  the  solar  and  the 
sensor’s  zenith  angles  while  the  relative  azimuth  an¬ 
gle  (^)  is  kept  constant  at  90®,  atj  which  in-sit u  above¬ 
water  measurements  are  usually  performed  in  order 
to  minimize  sky  and  Sun  glint.  i  In  the  figure,  varia¬ 
tions  in  both  contour  structure  and  color  intensity 
are  observed  for  the  coefficient  images  of  the 
different  wavelengths,  thereby  Confirming  the  wave¬ 
length  dependency  of  the  Rrs-<i;  relationship. 
Figure  4  also  shows  the  variability  of  the  coefficients 
a,  with  respect  to  <pix  axis)  and  0g  (y  eixis),  where  is 
kept  constant  at  40®  for  the  same  wavelengths  as 
in  Fig.  3. 

4.  Assessment  of  the  Bidirectional  Correction 
Algorithm 

4.1 .  Theoretical  Analysis 

Based  on  the  simulated  dataset,  we  carried  out  the 
assessment  of  the  bidirectional  Correction  performed 
with  the  proposed  model  (hereafter  referred  to  as 
the  CCNY  algorithm)  as  well  as  with  the  current 
operational  algorithm  [51  (referred  to  as  the  MG 
algorithm).  It  should  be  noted  that  although  the 
MG  algorithm  is  used  operationally  by  current 
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Fig.  2.  (Color  online)  Rr8(A)  with  respect  to  o)(>l)  at  viewing  angle  9^,  =  40°  for  X  =  551  nm  (top  row)  and  X  =  41^  nm  (bottom  row)  and  for 
the  relative  azimuth  angles  4>  =  45°,  90°,  and  180°  (from  left  to  right).  Simulations  are  shown  for  the  three  solar  aienith  angles  9^  =  0°  (red), 
30°  (green),  and  60°  (blue).  I 


above-water  and  satellite  data  processing  [40]  for  bi¬ 
directional  effect  correction  regardless  of  the  water 
type,  it  was  initially  developed  for  case  1  water  con¬ 
ditions.  Thus,  this  theoretical  analysis  mainly  in¬ 
tends  to  identify  the  amount  of  deviation  induced 
by  the  case  1  optimized  bidirectional  algorithm  in  ty¬ 
pical  coastal  water  conditions.  In  the  following  ana¬ 
lysis  we  focus  on  Rrs(>l)  for  a  fixed  viewing  angle 
{0^  =  40®)  at  which  the  sea-based  above-water  radi¬ 
ance  measurements  are  usually  performed  [41,42]. 
The  solar  zenith  angle  ranges  from  20°  to  70°,  and 
the  relative  azimuth  angle,  from  60°  to  180°.  These 
solar  zenith  and  relative  azimuth  angle  ranges  cover 
the  realistic  conditions  in  which  above-water  mea¬ 
surements  can  be  performed  with  the  relatively 
low  contamination  from  environmental  effects  [151. 

In  our  assessments,  simulated  remote  sensing  re¬ 
flectance  at  various  viewing  and  illumination  geome¬ 
tries,  Rrs(6?i,,  Og,  (f),  A)y  is  initially  translated  into  nadir 
remote  sensing  reflectance,  Rrs(0y  =  0, =  0, 
(^  =  0,  A),  using  lx)th  CCNY and  MG  algorithms.  Then 
those  nadir  remote  sensing  reflectance  values. 


Tabki  2.  Coefficients  a/,  sr^ ,  for  the  Nadir  Viewing  and  Solar  Angles 


A(mn) 


412 

443 

491 

551 

668 

«1 

«2 

«3 

0.04115 

0.10186 

0.01986 

0.04073 

0.09755 

0.002898 

0.04030 

0.11102 

0.00524 

0.04023 

0.1171 

-0.01355 

0.03879 

0.09425 

0.09660 

=  0°  and  ^5  =  0°. 


Rrs(0y  =  0,  =  0,  (^  =  0,  A),  obtained  with  the  CCNY 

algorithm  (from  here  on  denote(i  as 
and  the  MG  algorithm  (denoted  as  RrSretrieved(MGr 
compared  to  the  exactly  corresponding  simSated  na¬ 
dir  remote  sensing  reflectancO  value,  denoted  as 
RrsLuai-  Rrs?etrieved(CCNY)  Calculated  as  follows: 


0,035  0,04  0,045 


aj(443nm) 


a^(443nm) 


■0.2  -O.l  0 
a^(443nmj 


0  20  40  60  80  0  20  40  60  80  0  20  40  60  80 

(degree)  0^.  (degree)  0^  (degree) 


80 

?  60 

I  40 
<b'  20 


Fig.  3.  (Color  online)  Color  plots  of  the  coefficients  ai(9g ,  A), 

sr’^:  ai  (first  column),  02  (second  column),  and  (third  column) 
are  shown  as  a  function  of  the  sensoF^  zenith  angle  9^  {x  axis) 
and  the  solar  zenith  angle  (y  axis)  While  the  relative  azimuth 
angle  <t>  is  kept  constant  at  90°.  The  Coefficients  are  shown  for 
two  wavelengths:  X  =  443  nm  (top  row)  and  551  nm  (bottom  row). 
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(i)  a){X)  is  derived  by  fitting  to  the 

respective  «,{«?„  6)„,  <p,  X)  of  Eq.  (4);  (ii)  Rrs^j^evedcccNY) 
is  then  calculated  by  plugging  obtained  in  the 
first  step  into  Eq.  (4)  again,  together  with  =  0, 
Ov  =  0,<f)  =  0,A). 

However,  in  the  case  of  the  MG  algorithm,  prior  es¬ 
timation  of  [Chi]  is  required.  This  estimation  of  [Chi] 
is  carried  out  by  using  an  empirical  formula  that  re¬ 
lates  the  ratio  of  the  water-leaving  radiances  of  the 
blue  and  green  bands  to  [Chi]  [5,^] .  In  the  MG  algo¬ 
rithm  considered  in  this  study,  we  used  the  empirical 
formula  (OC2)  developed  as  a  product  of  the  Sea- 
WiFS  Bio-optical  Algorithm  Mini  Workshop  (Sea- 
BAM)  [43],  This  empirical  formula  is  also  adopted 
in  the  operational  AERONET  Ocean  Color  (OC)  data 
processing  for  ocean  color  satellite  data  validation 
activities.  The  formula  is  defined  as 


where 

(18) 

X  =  logic  1 

mrs(490)\ 

Urs(555)j' 

(19) 

The  use 

of  this 

blue-green  ratio 

[Chi] 

estimation  algorithm  [Eq.  ( 1B)1  remains  highly  ques¬ 
tionable  for  the  optical  properties  typically  retrieved 
in  coastal  waters  [44,45].  Nevertheless,  the  focus  of 
our  analysis  is  to  identify  the  deviation  induced  by 
the  whole  operational  bidirectional  correction  pro¬ 
cess.  Thus,  this  algorithm  has  been  used  here  with¬ 
out  modifications  or  adjustments.  On  the  other  hand, 
even  if  [Chi]  values  are  retrieved  with  reasonable  ac¬ 
curacy,  the  current  operational  MG  algorithm  will 
still  suffer  from  its  single  particle  type  assumption, 
in  which  the  contributions  of  scattering  from  the 
mineral  rich  NAP  (which  are  typical  for  coastal  water 


a|(443nm)  a2(443nm)  a^(443nm) 


a,(55lnm>  0^(55  Inm)  aj(55liim> 


0  60  120  180  0  60  120  180  0  60  120  180 


♦(degree)  ♦(degree)  ♦(degree) 

Fig.  4.  (Color  online)  ai  (first  column),  02  (second  column),  and  ^3 
(third  column),  sr"',  are  shown  as  a  function  of  4>  (x  axis)  and  the 
solar  zenith  angle  6,  iy  axis)  while  the  sensor’s  viewing  angle  is 
kept  constant  at  40^.  Color  intensity  values  are  the  same  as  in 
Fig.  3. 


conditions)  are  not  considered,  lix  our  simulations  we 
assumed  for  Eq.  (19)  Rrs(490) » Rrs(491)  and 
Rrs(555)  w  Rrs(551),  which  almost  did  not  affect 
[Chi]  values. 

Based  on  the  estimated  [Chi]  [Eq.  (18)],/'/Q  factors 
are  acquired  from  the  lookup  t^bles"T5l.  After  that, 
RrSretrieved(MG)  Calculated  as 


Rrs® 

^  “retrieved  (MG) 


JRo 


5R(0„,W) 


./•o(/t,Chl) 

‘Qo(>l,Chl) 

(^,/l,Chl) 


A0„^,Chl) 


xRrs{d„0u,^,X),  (20) 


where  SRq  denotes  the  particular  value  of  JR  when 

=  0.  In  the  same  way,  fo  denotes  the  particular 
value  of  f  when  6,  =  0,  and  Qo  denotes  the  Q  value 
when  <?(,  =  0  and  6,  =  0.  It  is  wjarth  noting  that  the 
parameters  f  and  Q  of  Eq.  (^)  are  computed  based 
on  typical  oceanic  case  1  water  optical  properties. 

Then,  both  Rrs^trieved(CCNY)  and  Rrs^^^ev^cMG)  are 
compared  to  Rrs^^^^,  and  the  following  statistical  in¬ 
dicators  are  calculated.  The  absolute  percent  differ¬ 
ence  for  each  matchup,  denoted  as  5,,  is 


5i  =  100  X  .  (21) 

Xi 

The  absolute  average  percent  difference  (AAPD)  is 
defined  as 


^D  =  iE‘5.’  (22) 

1  =  1 

and  the  average  percent  difference  (APD)  is 

=  (23) 

where  x  stands  for  Rrs^^y^j  and  y  for  either 

Rrs^trieved(CCNY)  RrSretrieved(MG)-  I"  the  compari- 
sons,  the  statistical  indicator  AAPD  provides  the  in¬ 
formation  regarding  the  dispersion,  while  APD  can 
be  used  to  assess  the  expected  bias  between  the  com¬ 
pared  datasets. 

Those  comparisons  have  been  carried  out  for 
(i)  Rrs^^y^j  and  corresponding  Rrs  (before  the  bidir¬ 
ectional  effect  is  removed)  in  order  to  make  the  as¬ 
sessment  of  the  bidirectional  variations  of  the 
assumed  water  types  for  the  viewing  and  illumina¬ 
tion  configurations  mentioned  above  (i.e., 

20®^,  ^  70®,  and  60®$<^  ^  180®),  (ii)  and 

^^retrieved(MG)  ^  identify  the  relative  error  resulting 
from  the  use  of  the  case  1  water  optimized  correct¬ 
ion  algorithm  for  case  2  water  conditions,  and 
(iii)  RrsO^^^  and  Rrs^t^evedcccNY)  obtain  the  statis¬ 
tics  for  the  error  induced  in  the  correction  processes. 
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Figures  5(aMe)  show  the  comparisons  between 

^•^retrieved  ^actual  443,  491,  551,  and 

668  nm,  respectively,  for  the  viewing  and  illumina¬ 
tion  configurations  mentioned  above.  It  is  observed 
that  the  comparisons  show  high  AAPD  values 
(6.26%  to  9.5%)  and  APD  values  (2.8%  to  9.4%)  at 
all  the  wavelengths  with  the  MG  algorithm,  and  re¬ 
gression  line  slopes  between  and 

Rrs^ctuai  l^OSy  0.98,  0.93,  and  0,97,  respec¬ 

tively,  for  412,  443,  491,  551,  and  668  nm  wave¬ 


lengths.  The  dispersion  between  the  two  datasets 
is  observed  to  be  relatively  smaller  for  the  shorter 
wavelengths  (412  and  443  nm)  and  is  largest 
at  551  nm. 

By  contrast,  the  comparisons  between 
RrSretrieved(CCNY)  ^nd  Rrs“rtuai  exhibit  much  smaller 
AAPD  (0.429M).8%)  and  APD  (0.054^.55%)  at 
all  wavelengths.  The  regres$ion  lines  between 
RrSretrieved(CCNY)  ^nd  ate  very  close  to  the 

1 !  1  line  for  all  wavelengths.  Moreover,  the  retrieved 
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Fig.  5.  (Color  online)  Comparisons  between  RrsJetuaj  derived  with  MG  (blue  dots)  and  CCNY  (green  dots)  algorithms  for 

0^  =  40®,  20®  S  5  70®,  and  60®  S  ^  S  180®  at  (a)  412,  (b)  443,  (c)  491,  (d)  551,  and  (e)  668  nm.  Regression  lines  between  Rrs^^rievedtMC)  ^d 
R™2ctu.i  are  shown  in  red,  and  those  of  are  shown  in  black,  (f)  Distribution  of  the  absolute  percent  difference,  6,  values 

between  RrsJrtuAi  ^d  Rrs  without  BRDF  correction  (red),  (black),  and  (blue)  for  each  matchup. 
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differences  between  Rrs2,jrieved(CCNY)  and  Rrs^^^, 
which  mostly  originate  from  the  least-mean-square 
fitting  procedure,  are  found  to  be  negligible  in  com¬ 
parison  with  the  actual  bidirectional  variation, 
which  is  generally  in  the  range  of  10%-20%,  as  can 
be  readily  seen  in  Fig.  5(f).  In  this  figure,  the  histo¬ 
gram  of  the  absolute  percent  difference  values 
combined  for  all  wavelengths  shows  that  d  values  be¬ 
tween  Rrs^^^^j  and  Rrs  without  BRDF  correction  are 
in  the  range  of  2.3%-25%  and  AAPD  is  around 
13.5%,  which  is  consistent  with  previous  studies 
[6,7].  In  addition,  it  can  be  seen  that  without  the 
BRDF  correction,  only  9.5%  of  the  data  have  absolute 
percent  difference  (6)  values  within  5%,  and  more 
than  two-thirds  of  the  data  have  <5  values  greater 
than  10%.  These  results  demonstrate  the  importance 
of  the  bidirectional  effect  correction  in  the  retrieval 
process  if  the  water-leaving  radiances  (and  remote 
sensing  reflectances)  are  to  be  retrieved  within  the 
ocean  color  sensor  community’s  targeted  accuracy  le¬ 
vel  (i.e.,  less  than  5%  uncertainty)  for  above-water 
measurements. 

Similarly,  S  values  between  Rrs®^^^  and 
^^®retrieved(MG)  from  0%  to  27%.  Nevertheless, 

AAPD  is  reduced  to  7.56%.  Moreover,  it  was  observed 
that  43%  and  70.6%  of  the  data  set  have  S  values  less 
than  5%  and  10%,  respectively.  These  results  show 
that  the  bidirectional  effect  is  reduced  with  the 
use  of  the  MG  correction  algorithm  and  that  using 
the  MG  algorithm  in  coastal  water  conditions  can 
still  offer  better  retrievals  than  not  correcting  at 
all.  However,  57%  of  the  data  have  <5  values  greater 
than  5%,  and  thus  these  results  also  outline  the  need 
for  the  bidirectional  correction  algorithm  optimized 
for  case  2  water  conditions.  On  the  other  hand,  the 
distribution  curve  of  S  values  between  Rrs^^^y^j  and 
^®retrieved(CCNY)  shows  a  much  narrower  range, 
and  98.4%  of  the  data  have  <5  values  less  than  5%. 
The  mean  of  the  <5  value  is  also  drastically  reduced 
to  0.69%,  demonstrating  that  the  CCNY  BRDF  cor¬ 
rection  algorithm  works  very  well  with  the  simulated 
dataset. 

4.2.  Noise  Sensitivity  Analysis 

Water-leaving  radiances  retrieved  from  above-water 
radiometric  measurements  are  subject  to  uncer¬ 
tainty  that  arises  from  geophysical  noise,  such  as  ra¬ 
pid  changes  in  environmental  conditions,  and  Sun 
and  sky  glints,  which  are  randomly  fluctuating  from 
the  effects  of  surface  waves.  The  current  above-water 
radiometry  processing  algorithm  employs  the  glint 
filtering  method,  in  which  only  the  lowest  of  multiple 
successive  total  sea  radiance  measurements  taken 
within  a  short  period  are  utilized  in  the  processing 
of  Rrs  in  order  to  minimize  the  effect  of  the  glint 
[1^,41].  Such  practices  may  result  in  the  residual 
random  component  affecting  the  retrieved  Rrs.  Re¬ 
moval  of  the  reflected  sky  radiance  can  also  lead 
to  residual  errors,  which  are  not  taken  into  account 
by  the  removal  model  [1^].  For  the  satellite  data,  er¬ 
rors  can  arise  from  inaccuracies  in  atmospheric  cor¬ 


rection  or  from  other  sources.  The  analysis  of  such 
residual  errors  resulting  from  each  step  of  the  data 
processing  procedures  is  beyond  the  scope  of  this  pa¬ 
per  and  deserves  a  dedicated  study  because  of  the  im¬ 
portance  of  the  topic.  Here,  Our  analysis  simply 
intends  to  identify  the  proposed  model’s  sensitivity 
to  the  noise,  in  other  words,  to  verify  whether  the 
proposed  model  can  still  offer  meaningful  corrections 
(i.e.,  whether  output  Rrs^tneved  closer  to  corre¬ 
sponding  Rrs^^^^j  than  Rrs(0t^,0,,<^,  A)  are)  when 
Rrs  measurements  are  contaminated  by  geophysical 
noise,  and  to  determine  the  threshold  of  noise  level 
with  respect  to  the  Rrs  signal.  Since  the  proposed 
model  performs  the  correction  at  each  wavelen^h  in¬ 
dependently  and  modeling  error$  at  each  wavelength 
are  much  lower  than  assumed  noise  levels,  the  mod¬ 
el’s  sensitivity  to  the  noise  at  each  wavelength  can  be 
evaluated  from  this  analysis. 

In  order  to  test  the  applicability  of  the  CCNY 
algorithm  in  the  above-water  correction  procedures, 
we  first  added  synthetic  noise  signals  to  the 
values  before  the  BRDF  correction 
was  applied.  Then,  these  noise-contaminated 
Rrs{0jj,0s,(p,X)  values  are  converted  to  nadir  remote 
sensing  reflectances,  using  the  CCNY  algorithm. 
After  that,  noise-contaminated  Rrs^gy^^^CNY) 
lues  are  compared  to  the  respective  Rrs^^^^  values 
as  in  the  previous  section.  For  this  analysis,  s)ntithetic 
noise  signals  are  independently  generated  for  each 
Rrs(^y,0,,<^,  A)  value  as  random  numbers  with  a 
Gaussian  distribution  with  zero  mean  and  standard 
deviation  proportional  to  the  corresponding  Rrs 
values. 

Figure  6  shows  the  distributions  of  5  values  (com¬ 
bined  for  all  wavelengths)  between  Rrs9 


noise-contaminated  RrSretrieYed( 
levels  of  noise  (5%,  10%,  and  1 


^actual 

values  for  three 
yc)  together  with  the 
distribution  of  <5  values  for  Rrs  without  a  BRDF  cor¬ 
rection.  It  was  found  that  AAPD  values  of  the  retrie¬ 
vals  are  increased  to  4.0%,  7.9%,  and  12.9%, 
respectively,  for  the  dataset  with  5%,  10%,  and 
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Fig.  6.  (Color  online)  Distribution  of  the  absolute  percent  differ¬ 
ence  (S)  values  for  inff&ctcd  with  four  different 

noise  levels:  0%  (blue),  5%  (green),  10%  (black),  and  15%  (magen¬ 
ta).  The  red  curve  corresponds  to  the  cAse  where  Rrs  is  not  cor¬ 
rected  for  BRDF.  Illumination  and  vi0wing  configurations  are 
the  same  as  in  Fig.  5. 
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15%  noise  levels.  Although  retrieval  performances 
are  lower  with  the  introduction  of  noise  to  the  signal, 
<5  values  for  all  noise  levels  are  shifted  toward  smaller 
values  compared  to  the  S  values  for  Rrs  without 
BRDF  correction.  Thus,  it  shows  that  even  noise- 
contaminated  Rrs  measurements  can  still  consider¬ 
ably  benefit  from  bidirectional  correction. 

5.  Bidirectional  Reflectance  Distribution  Correction 
of  Above-Water  Measurements 

5.1.  LISCO  instrumentation 

As  was  mentioned  above  in  the  introduction,  perfor¬ 
mance  of  the  proposed  BRDF  correction  approach 
was  also  evaluated  by  measurements  obtained  at 
LISCO.  LISCO  combines  a  multispectral  SeaPRISM 
(Cimel,  France)  instrument  as  a  part  of  AERONET- 
OC  Network  [461  with  a  colocated  hyperspectral 
HyperSAS  (Satlantic,  Canada)  radiometric  system. 
Both  instruments  are  positioned  on  a  retractable 
tower  with  an  elevation  of  12  m.  They  were  installed 


are  achieved  for  180  channels  regularly  spaced 
between  305  and  905  nm.  Tvfo  radiance  sensors 
are  moimted  on  the  same  instalment  plate  as  the 
SeaPRISM  system  and  have  a  fiill  angle  field  of  view 
of  3®.  But,  unlike  SeaPRISM,  the  sensors  point  ex¬ 
actly  westward  all  the  time.  As  a  result,  the  relative 
azimuth  angle  0  changes  throughout  the  day,  encom¬ 
passing  the  whole  range  from  0^  to  180®.  HyperSAS 
measurement  sequences  are  also  executed  every 
30  min,  as  in  the  case  of  SeaPRISM,  and  each  se¬ 
quence  takes  2  min.  Within  2  min,  the  number  of 
measurements  acquired  by  H)q)arSAS  system  varies 
from  45  to  80  for  and  from  50  to  220  for  Lg^y  and 
measurements  depending  on  ambient  light 
conditions. 

5.2.  Method  and  Data  Processing 

Remote  sensing  reflectance  spectra,  Rrs(A),  for  a  spe¬ 
cific  viewing  and  illumination  irondition  are  calcu¬ 
lated  as  described  in  [16]  as  follows: 


Rrs(^g,^?t;,0,A) 


Ea{e,,X) 


(24) 


on  the  LISCO  platform  in  October  2009  and  have 
been  continuously  providing  data  since  then. 

The  SeaPRISM  system  consists  of  CE-318  Sun 
photometers  modified  to  meet  requirements  for 
above-water  radiometry.  The  photometers  perform 
radiance  measurements  with  a  full-angle  field  of 
view  of  1.2®  to  determine  the  total  radiance  from 
the  sea,  and  the  sky,  A), 

where  =  180  -  0^.  Thanks  to  the  rotating  feature 
of  SeaiPRISM,  0  is  always  set  to  90®  regardless  of 
the  Sun  position,  and  O^j  is  set  at  40®.  These  values 
were  determined  in  order  to  minimize  the  perturba¬ 
tions  resulting  from  the  Sun  glint  off  the  sea  surface 
[16]  and  the  deployment  superstructure  itself  or  its 
shadow  [411.  The  SeaPRISM  system  configuration  of 
LISCO  performs  ocean  color  measurements  at  the 
412,  443,  491,  551,  668,  870,  and  1018  nm  center  wa¬ 
velengths.  These  center  wavelengths  were  selected  to 
be  as  close  as  possible  to  the  bands  of  current  ocegm 
color  earth  observation  systems  in  order  to  support 
essential  validation  activities.  Measurement  se¬ 
quences  are  executed  every  30  min  within  it4h  of 
12:00  PM  local  time.  Each  measurement  sequence, 
which  consists  of  three  sky  (Lgi^y)  and  eleven  sea 
{Lf)  radiance  data  acquisitions,  takes  approximately 
6  min. 

The  hyperspectral  measurements  are  achieved  by 
a  H3perSAS  system  providing  high  precision  mea¬ 
surements  of  L8i^(^g,^J,0,A),  and 

down  welling  spectral  irradiance  Ed{0^,X),  The  radi¬ 
ance  and  irradiance  measurements  of  H)perSAS 


Following  our  previous  study  [15],  the  total  radiance 
(Lj)  utilized  in  Eq.  (24)  for  the  SeaPRISM  instrument 
is  obtained  from  Level  1.5  cloud  screened  data  of 
AERONET-OC  data  distribution  [^1.  In  the  case 
of  H)perSAS,  is  calculated  by  averaging  the  lowest 
5%  of  the  sea  radiance  measurepients  taken  during 
the  2  min  data  acquisition  time.  H5perSAS  performs 
45  to  80  total  sea  radiance  measurements  during  its 
2  min  data  acquisition  window,  depending  on  the 
brightness  of  the  total  radiances  being  obtained.  Tak¬ 
ing  the  lowest  5%  allows  the  averaging  of  at  least  two 
measurements  during  every  cyde  of  2  min  data  ac¬ 
quisition  for  further  processing.  On  the  other  hand, 
the  sky  radiance,  Lg,  is  determined  by  simply  aver¬ 
aging  all  sky  radiance  measurements.  The  sea  sur¬ 
face  reflection  coefficient,  p{d^,0^j^<p,W)y  for  the 
desired  exact  viewing  and  illumination  geometries 
and  the  wind  speed  are  obtained  by  performing  a  lin¬ 
ear  interpolation  on  the  reflection  coefficient  dataset 
acquired  through  Hydrolight’s  simulations  described 
in  Section  3.  Input  surface  wind  speed  data  for  each 
data  acquisition  cycle  are  acquired  from  the  National 
Weather  Service.  Downwelling  irradiance,  Erf(A), 
measurements  are  acquired  by  the  vertically  in¬ 
stalled  irradiance  sensor,  part  of  the  H3perSAS  sys¬ 
tem  and  used  in  the  calculation  of  Rrs  (A)  for  both 
H3perSAS  and  SeaPRISM  systems. 

5.3.  Assessment  of  the  Bidirectional  Reflectance 
Distribution  Function  via  LISCO  Data 

One  of  the  major  advantages  of  the  LISCO  instru¬ 
mentation  set  is  that  HyperSAS  and  SeaPRISM 
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instruments  acquire  data  almost  concurrently  but 
with  different  viewing  geometries  as  the  day  pro¬ 
gresses,  except  for  the  time  when  the  Sun  is  exactly 
south  when  both  instruments  point  west.  As  both  in¬ 
struments  make  measurements  within  a  10  min  win¬ 
dow,  it  is  reasonable  to  consider  that  the  composition 
of  the  water  constituents  within  that  time  window 
remains  the  same.  Ideally,  observed  variations  in 
the  water-leaving  radiances  acquired  by  the  two  in¬ 
struments  can  be  attributed  to  the  variations  due  to 
the  BRDF  only.  Hence,  the  assessments  of  the  perfor¬ 
mance  of  the  BRDF  corrections  can  be  made  by  com¬ 
paring  the  nadir  remote  sensing  reflectance  obtained 
from  SeaPRISM  [denoted  as  RrSspj.(A)]  and  from  Hy- 
perSAS  [denoted  as  Rrs^(A)].  However,  it  should  1^ 
pointed  out  that  directional  fluctuations  in  the  mea¬ 
sured  sea  radiance  are  also  induced  by  the  Sun  and 
sky  light  reflections  on  the  ruffled  surface,  which  in 
the  end  may  create  significant  uncertainties  in  the 
retrieved  water-leaving  radiance  [15,48].  Therefore, 
it  is  of  paramount  importance  to  filter  out  the  data 
corrupted  by  unexpected  environmental  effects  or 
any  stochastic  artifacts  beforehand.  The  filtering 
procedure  for  quality  assurance  has  already  been  de¬ 
veloped  for  the  SeaPRISM  system  in  the  framework 
of  AERONET-OC  distribution  [46],  and  all  the  data 
we  used  in  this  study  are  quality-checked  level  1.5 
cloud-screened  data  [47].  Similarly,  a  specific  Hyper- 
SAS  data  filtering  procedure  has  been  developed  by 
us  independent  of  SeaPRISM  [15],  and  those  proce¬ 
dures  were  closely  followed  in  this  study  for  the  data 
processing.  Then,  hyperspectral  HyperSAS  data 
were  integrated  with  the  sensor  relative  spectral  re¬ 
sponse  fimction  of  each  SeaPRISM  band  in  order  to 
produce  equivalent  data  for  both  systems.  The  data 
involved  were  restricted  to  SeaPRISM  measurement 
sequences  taken  within  ±10  min  of  HyperSAS  se¬ 
quence  intervals.  In  addition,  findings  from  our 
previous  study  suggest  that  the  contamination  by  en¬ 


vironmental  effects  of  the  retrieved  water-leaving  ra¬ 
diance  becomes  sensitive  for  (p  <  60®  regardless  of 
the  Sun^s  elevation,  with  uncertainties  higher  than 
5%  [1^].  As  a  consequence,  it  has  been  decided  to 
eliminate  HyperSAS  measurements  taken  with  p  < 
60®  from  the  analysis. 

In  the  following,  we  first  show  the  intercompari¬ 
sons  between  the  remote  sensing  reflectance  of 
HyperSAS  and  SeaPRISM  before  the  bidirectional 
effect  are  corrected  [denoted  as  RrsHs(^)  and 
RrssprW>  respectively].  This  coipparison  (see  Fig.  7) 
allows  us  to  make  the  assessment  of  the  variation  in 
two  measurements  before  the  bidirectional  effects 
are  removed.  Second,  we  perform  the  comparisons 
between  the  final  nadir  remote  sensing  reflectance 
of  both  instruments  [denoted  as  Rrs^g(A)  and 
RrsgprfA)]  processed  with  both  the  CCNY  and  MG  al¬ 
gorithms  (see  Fig.  8).  Finally,  we  compare  the  final 
nadir  remote  sensing  reflectance  of  both  instruments 
processed  with  the  CCNY  algorithm  (denoted  as 
^^retrieved(CCNY)^  ^  those  processed  with  the  MG  al- 
gorithm  (denoted  as  Rrs^trievedono))  2)- 

the  comparisons,  becaxise  neither  SeaPRISM  nor  Hy¬ 
perSAS  data  can  be  assigned  as  the  reference,  the  fol¬ 
lowing  statistical  indicators  given  in  percentage 
points  are  used. 

The  absolute  relative  percent  difference  (ARPD)  is 
defined  as  follows: 


ARPD  =  2^flaz£il, 


(25) 


and  the  unbiased  relative  percent  difference  (URPD) 
is 


URPD  =  ^ 

N  ^Xi+yi 


(26) 


•  4I2nni  •443nni  491  nm  •35lnm  •  668nni 

(a)  (b) 


Fig.  7.  (Color  online)  Intercomparisons  of  remote  sensing  reflectances  (in  sr“' )  derived  from  SeaPRISM  and  HyperSAS  before  correction 
for  the  bidirectional  effect:  (a)  relative  azimuth  angles  for  HyperSAS  observations  are  restricted  in  the  60®  ^  180®  range;  (b)  relative 

azimuth  angle  range  is  restricted  to  80®  100®.  N  is  the  total  number  of  the  comparisons,  and  the  value  in  parentheses  is  the  number  of 

different  measurement  sequences  used  in  the  comparison.  | 
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Fig.  8.  (Color  online)  Intercomparisons  of  SeaPRISM  and  HyperSAS  remote  sensing  reflectance  measurements  (in  sr~')  after  the  bidir¬ 
ectional  effect  is  corrected:  (a)  processed  with  the  MG  algorithm,  (b)  processed  with  the  CCNY  algorithm.  Relative  azimuth  angles,  0,  for 
HyperSAS  observations  are  the  same  as  in  Fig.  7(a>. 


where  SeaPRISM  data  are  taken  as  x  and  HyperSAS 
data  are  taken  as  y. 

Figure  7(a)  shows  the  comparison  of  remote  sen¬ 
sing  reflectance  measurements  of  HyperSAS  and 
SeaPRISM  [RrsHsC'^)  and  Rrsspr('i)]  before  the  bidir¬ 
ectional  effect  is  removed  for  afl  the  data  used  in  this 
analysis.  The  comparison  exhibits  strong  correla¬ 
tions  between  HyperSAS  and  SeaPRISM  data,  with 
the  value  around  0.967,  and  the  regression  line 
between  two  pieces  of  data  is  also  very  close  to  the 
1 : 1  line.  The  URPD  and  ARPD  values,  which  show 
the  bias  and  dispersion  between  two  pieces  of  data, 
are  8.73%  and  12.3%,  respectively.  However,  when 
the  relative  azimuth  angle  range  for  the  HyperSAS 
instrument  is  restricted  to  the  80°  ^  100°  range 

[see  Fig.  7(h)l,  the  URPD  and  ARPD  values  are  re¬ 
duced  down  to  -0.342%  and  6.02%,  respectively  In 
addition,  the  value  is  increased  to  0.991,  and 
the  regression  line  comes  closer  to  the  1 1 1  line,  show¬ 
ing  a  better  consistency  between  the  two  measure¬ 
ments.  In  fact,  with  this  restricted  relative 
azimuth  range,  HyperSAS  and  SeaPRISM  measure¬ 
ments  are  taken  for  an  azimuth  within  ±10°,  and 
thus  bidirectional  effects  for  the  measurements  of 
both  systems  are  about  the  same.  Therefore,  it  can 
be  considered  that  this  statistical  difference  between 
the  two  comparisons  arises  from  the  bidirectional 
variation  of  the  water-leaving  radiance  as  well  as 
from  possible  residual  Sun  and  sky  glint  effects.  In 
fact,  these  residual  glints  may  be  considered  to  have 
the  same  effect  as  the  noise  on  the  correction  proce¬ 
dures,  as  shown  in  Section  4.2. 

Comparisons  of  the  Rrs£^)  and  RrSgpj.(A)  pro¬ 
cessed  with  the  MG  and  CCNi  algorithms  are  shown 
in  Figs.  8(a)  and  8(b),  respectively.  In  the  comparison 
of  the  Rrs^(A)  and  Rrs^pj.(A)  processed  with  the  MG 
algorithm,  the  URPD  value  is  reduced  down  to 
7.67%,  showing  a  1.06%  reduction  in  relative  percent 
difference  compared  to  the  same  data  not  corrected 
for  the  bidirectional  effects.  However,  its  ARPD  value 
is  increased  by  0.3%  compared  to  the  value  before  the 


correction.  The  B^  value  is  also  observed  to  be  re¬ 
duced  down  to  0.958.  These  results  indicate  that 
the  dispersion  between  the  two  data  is  somewhat  in¬ 
creased  by  the  MG  correction  procedure.  This  shows 
the  inadequacy  of  bidirectional  correction  based  on 
the  assumption  of  case  1  water  for  the  t3npical  LISCO 
waters.  On  the  other  hand,  comparison  with  the 
same  data  processed  with  the  CCNY  algorithm 
shows  a  reduction  in  both  ARPt)  and  URPD  values 
(2.1%  and  3.14%,  respectively)  with  respect  to  the 
data  not  corrected  for  bidirectional  effects,  while 
B^  is  increased  to  0.972.  Therefore,  this  indicates 
that  with  the  use  of  the  CCNY  bidirectional  correc¬ 
tion  algorithm,  dispersion  between  the  two  mea¬ 
surements  resulting  from  the  different  viewing  con¬ 
figuration  can  be  reduced,  thereby  demonstrating 
the  relevance  and  appropriateness  of  the  CCNY  bi¬ 
directional  correction  for  the  water-leaving  radiance 
measured  by  the  H3nperSAS  and  SeaPRISM  instru¬ 
ments  in  the  moderately  turbid  waters  of  the  LISCO 
site.  It  can  also  be  observed  that  both  ARPD  and 
URPD  values  of  the  data  corrected  with  the  CCNY 
algorithm  are  lower  than  those  corrected  with  the 
MG  algorithm  at  every  wavelength,  showing  1.95% 
to  3.0%  improvement  in  terms  of  ARPD  and  1.65% 
to  2.65%  in  terms  of  URPD  (see  Table  3  for  details). 

Figure  9  shows  the  intercomparisons  between 
R'^retrievedcccr^  and  Rrsl^trie^^(uGy  It  is  observed 
that  the  relative  absolute  percent  difference  between 
the  output  nadir  remote  sensing  reflectance  of  the 
two  algorithms  ranges  from  0%  to  14.5%,  where 
the  largest  difference  found  at  551  nm  with  the 
AAPD  between  the  outputs  of  the  two  algorithms 
is  6.14%.  The  slopes  between  the  output  final  nadir 
remote  sensing  reflectances  of  the  two  algorithms  are 
1.05,  0.95,  and  0.905  for  412,  49l,  and  551  nm  wave¬ 
lengths,  respectively,  suggesting  that  the  MG  algo¬ 
rithm  is  overestimating  the  nadir  remote  sensing 
reflectances  at  412  nm  and  underestimating  them 
at  491  and  551  nm.  It  should  be  noted  that  compar¬ 
isons  carried  out  with  the  field  measurement  data 


232  APPLIED  OPTICS  /  Vol.  51 ,  No.  2  / 10  January  2012 


•  412  nm  •  443  nm 


491  nm 


•  55 1  nm  •  nm 


\  10*'  (a) 


^  lo  '  (b) 


Fig.  9.  (Color  online)  Intercomparisons  between  ^^retrieved(MG) 

and  (c)  668  nm.  (f)  Distribution  of  absolute  percent  difference  (^)  between  and  for  all  wavelengths.  Both 

HyperSAS  and  SeaPRISM  data  are  plotted  together. 


illustrate  similar  slopes  to  those  carried  out  with  the 
simulated  data  set. 

6.  Application  of  the  BRDF  Correction  Algorithm  to 
MODIS  Satellite  Data 

A  previous  study  has  confirmed  that  the  LISCO  site 
is  appropriate  for  use  in  calibrationA^alidation  of  the 
ocean  color  satellites  in  coastal  waters  [141.  In  the 
study,  we  compared  time  series  normalized  water¬ 
leaving  radiance  (nLw)  measurements  obtained  at 
LISCO  with  the  ones  derived  from  satellite  data. 
That  time  series  comparison  performed  with  the 
use  of  a  nine-month  data  period  exhibited  both  qua¬ 
litative  and  quantitative  agreements  between  water¬ 
leaving  radiances  derived  from  LISCO  measure¬ 
ments  and  satellite  data,  showing  strong  correlations 
and  relatively  low  average  percent  differences  be¬ 
tween  the  datasets,  especially  at  443,  491,  and 


551  nm  [14].  Similar  to  this  previous  study,  matchup 
comparisons  have  been  carried  out  between  the 
in  situ  SeaPRISM  and  satellite  Rrs^(A)  data  for  LIS- 
CO's  location.  Rrs^(A)  data  of  both  in  situ  and  satel¬ 
lite  types  are  processed  in  the  following  two  different 
ways:  (i)  with  standard  BRDF  correction  (the  MG  al¬ 
gorithm)  and  (ii)  with  the  proposed  CCNY  algorithm. 
Then,  in  situ  and  satellite  data  matchup  comparisons 
are  carried  out  for  both  obtained  datasets  in  order  to 
assess  merits  and  benefits  of  the  proposed  algorithm 
over  the  standard  one. 

We  used  the  data  from  the  MODIS  onboard  the 
Aqua  satellite  for  the  matchup  and  for  the  compari¬ 
sons  with  the  in  situ  data.  MODIS  (Aqua)  Level  1 
data  for  a  14-month  period  (October  2009  through 
December  2010)  were  acquired  from  the  ocean  color 
data  distribution  website  [49].  These  data  were  pro¬ 
cessed  using  the  software  package  SeaDAS  version 
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Fig.  10.  (Color  online)  Scattering  plots  of  the  comparisons  between  MODIS  and  in  situ  data;  (a)  compari$on  between  MODIS  and 
SeaPRISM  with  MG  BRDF  processing,  (b)  comparison  between  MODIS  and  SeaPRISM  with  CCNY  BRDF  processing. 


6.2  (which  contains  all  of  the  final  processing  pro¬ 
grams  and  the  source  code  used  for  the  OBPG’s 
Ocean  Color  Reprocessing  2010)  into  two  groups  of 
Level  2  data  files  (one  group  is  processed  using  the 
MG  BRDF  correction  and  another  group  with  the 
CCNY  BRDF  correction)  at  the  nominal  1  km  resolu¬ 
tion.  Then  we  extracted  the  data  of  the  3  km  x  3  km 
area  centering  the  LISCO  platform  to  create  two 
time  series  of  MODIS  Rrs^(A)  data  sets  processed 
with  the  MG  and  CCNY  BRDF  correction  algo¬ 
rithms.  For  data  quality  assurance  purposes,  Level 

2  flags  [i.e.,  land,  cloud,  high  Sun  glint,  very  low 
water-leaving  radiance  (cloud  shadow),  and  atmo¬ 
spheric  correction  failure]  were  applied  in  the  data 
extraction  procedure.  The  average  values  of  the 

3  km  X  3  km  area  are  considered  qualified  for  com¬ 
parison  with  in  situ  data  when  at  least  50%  of  the 
pixels  are  not  affected  by  these  standard  flags.  On 
the  other  hand,  two  sets  of  in  situ  SeaPRISM 
Rrs®(A)  data  processed  with  the  MG  and  CCNY 
BRDF  corrections  were  selected  from  the  measure¬ 
ments  that  passed  the  data  filtering  procedures 
[15]  and  taken  within  ±4  hours  of  the  satellite  over¬ 
pass  time  for  the  LISCO  location. 

Subsequently,  we  compared  the  MODIS  Rrs®(A) 
datasets  (processed  using  both  the  CCNY  and  MG  al¬ 
gorithms  for  bidirectional  correction),  from  here  on 


denoted  as  Rtsmodis^^)’  corresponding  SeaPR¬ 
ISM  (Rrsg  (A))  aatasets.  Comp^sons  were  carried 
out  for  both  correction  algorithm^  and  statistics,  such 
as  AAPD,  and  coefficients  of  determination  (R^)  have 
been  obtained.  A  total  of  42  MODIS  and  SeaPRISM 
matchups  and  comparisons  werO  carried  out.  Sensor 
and  solar  zenith  angles  of  MODIS  data  used  in  the 
comparisons  are  in  the  range  0f  5®  to  56®  and  20® 
to  60®,  respectively,  and  therefore  a  wide  range  of  il¬ 
lumination  and  viewing  geometries  are  available  to 
test  the  bidirectional  corrections.  Table  4  shows  the 
resulting  AAPD  values  from  the  comparison  between 
Rrs^ODis(^)  processed  with  the 

MG  ana  CCNY  algontnms. 

In  the  comparisons  of  Rrs^nis(^) 
data,  the  CCNY  algorithm  exmoite  noticeable  im¬ 
provements  throughout  the  spectrum.  The  spectrally 
averaged  AAPD  between  the  two  datasets  is  reduced 
about  3.15%  with  the  use  of  the  CCNY  algorithm  (see 
Table  4  for  details).  Figure  10  shows  the  scatterplots 
of  the  Rrs^oDis(^)  versusRrsQ_(jl)  data.  Figure  10(a) 
shows  the  comparison  of  MUDIS  and  SeaPRISM 
data,  processed  with  the  MG  algorithm,  while 
Fig.  10(b)  shows  the  comparison  of  the  same  data 
processed  with  the  CCNY  algorithm.  The  three 
wavelengths  491,  551,  and  668  nm  are  shown  in  the 
plots,  and  a  stronger  correlation  (0.926)  is  observed 


Tabia  3.  Statiaticai  Summary  of  the  intercompariaont  of  SeaPRISM  and  HyperSAS  Remote  Sensing  Reflectanca  Measurements 


Wavelength  (nm) 


412 

443 

491 

551 

668 

Spectral  Average 

Before  BRDF  correction 

R2 

0.7(0.81) 

0.85(0.94) 

0.95(0.98) 

0.96(0.98) 

0.95(0.05) 

0.967(0.991) 

ARPD 

21(11.3) 

16(5.81) 

7.86(3.23) 

5.9(3.12) 

10.8(6.^7) 

12.3(6.02) 

URPD 

12.1(-.5.9) 

11.6(-0  02) 

5.18(0.66) 

3.56(0.67) 

7.48(4. ’41) 

8.73(-0.34) 

Corrected  with  MG  algorithm 

0.689 

0.833 

0.941 

0.945 

0.883 

0.958 

ARPD 

20.35 

16.3 

8.17 

6.99 

11 

12.6 

URPD 

12.3 

11.8 

5.36 

3.76 

7.73 

7.67 

Corrected  with  CCNY  algorithm 

0.749 

0.878 

0.961 

0.957 

0.906 

0.972 

ARPD 

17.8 

13.3 

6.09 

5.04 

8.76 

10.2 

URPD 

9.65 

9.41 

3.46 

2.11 

5.28 

5.59 

First  row:  Before  the  bidirectional  effect  is  removed  (the  60®  S  ^  5  180®  range  is  given  first,  and  the  80®  ^4fS^  00®  range  is  shown  in  parentheses). 
Second  row:  After  the  bidirectional  effect  is  corrected  with  the  MG  algorithm.  Tlurd  row:  (Corrected  with  the  CCNY  algorithm. 
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Table  4.  AAPD  Valuea  from  the  Comparlaon  between  the 
RrsJiOQi^CA)  and  RrsJp,(A)  Data 

Wavelength  (nm) 


AAPD  (%) 

412 

443 

491 

551 

668 

MG 

46.43 

38.85 

16.68 

13.61 

24.54 

CCNY 

42.40 

34.16 

14.93 

10.99 

21.89 

Improvement 

4.03 

4.69 

1.75 

2.62 

2.65 

for  the  comparisons  processed  with  the  CCNY  algo¬ 
rithm.  The  dispersion  between  the  satellite  and  in 
situ  data,  esp)ecially  for  high  Rrs®(A)  values  (which 
are  usually  associated  with  high  particulate  loading 
in  the  water),  is  reduced  for  the  CCNY  processing, 
thereby  exhibiting  its  appropriateness  for  use  in  the 
bidirectional  correction  for  coastal  water  types. 

7.  Summary  and  Conclusion 

We  first  created  a  dataset  of  remote  sensing  reflec¬ 
tance  spectra,  Rrs(0„,0„  ^,A),  through  radiative 
transfer  simulations  made  with  Hydrolight.  Simula¬ 
tions  were  performed  for  all  ranges  of  viewing  and 
illumination  geometries,  with  input  parameters  of 
the  lOP  model  generated  as  uniformly  distributed 
random  variables  in  the  prescribed  ranges  typical 
for  coastal  waters.  Then,  Rrs(A)  values  simulated  for 
all  viewing  and  illumination  geometries  were  para¬ 
meterized  with  the  single  backscattering  albedo, 
a)(X)y  for  each  set  of  geometries.  Based  on  this  para¬ 
meterization,  we  proposed  an  LSQ  fit  third  order 
polynomial  remote  sensing  reflectance  model  as  a 
function  of  the  backscattering  albedo,  for  which  co¬ 
efficients  were  tabulated  as  a  lookup  table  for  the 
sets  of  illumination  and  viewing  geometries  and 
wavelengths  to  estimate  the  directional  variation 
of  the  water-leaving  radiance  field  of  t5npical  case  2 
water  conditions. 

We  then  carried  out  the  assessments  of  the  bidir¬ 
ectional  correction  performed  with  the  CCNY  model 
as  well  as  the  operational  MG  algorithm  based  on  the 
simulated  dataset.  Our  assessments  showed  that  the 
proposed  CCNY  algorithm  demonstrates  advantages 
over  the  operational  MG  algorithm.  The  retrieval  er¬ 
rors  resulting  from  CCNY  BRDF  corrections  are  re¬ 
latively  small  and  have  AAPD  values  in  the  range  of 
0.42%^.85%.These  small  retrieval  errors  mostly 
come  from  the  LSQ  fitting  procedure  and  are  well 
below  the  well-known  bidirectional  variation  of 
10%-20%.  In  addition,  the  regression  lines  between 
the  retrieved  and  the  actual  remote  sensing  reflec¬ 
tance  are  very  close  to  the  1 : 1  line  for  all  viewing 
and  illumination  geometries.  On  the  other  hand, 
comparisons  carried  out  using  the  current  MG  algo¬ 
rithm  depict  high  AAPD  values  in  the  range  of 
6.2%-9.54%.  In  addition,  more  than  57%  of  the 
MG  retrievals  are  outside  the  range  of  5%  accuracy 
targeted  by  the  OC  sensor  validation  procedure.  Sys¬ 
tematic  variations  of  the  slope  and  of  the  bias  with 
the  scattering  angles  are  also  observed  at  all  wave¬ 
lengths  for  the  comparisons.  These  results  indicate 


that  the  MG  algorithm,  which  is  optimized  for  the 
case  1  water  condition,  with  [Chi]  being  the  main  in¬ 
put,  is  unsuited  for  use  in  optimally  complex  case  2 
waters  if  the  OC  sensor  validation  procediu^^s  5%  re¬ 
trieval  accuracy  target  is  to  be  achieved. 

We  also  carried  out  comparisons  of  two  algorithms 
for  the  above-water  radiometric  dataset  of  LISCO, 
which  offers  the  unique  capability  of  making  near¬ 
concurrent  water-leaving  radiance  measurements 
from  different  viewing  geometries.  LISCO^s  instru¬ 
mentation  setup,  in  which  two  instruments  acquire 
the  data  almost  concurrently  but  at  different  scatter¬ 
ing  angles  with  respect  to  the  Sun,  is  the  major 
advantage  we  have  with  regard  to  analyzing  the  fi¬ 
delity  of  bidirectional  correction  algorithms:  we  can 
compare  the  final  remote  sensing  products,  Rrs®, 
retrieved  from  the  measurenxents  made  by  two 
instruments  acquiring  data  at  different  viewing  geo¬ 
metries.  Comparisons  are  performed  at  all  SeaPR- 
ISM  wavelengths.  With  the  proposed  algorithm, 
comparison  between  the  final  remote  sensing  reflec¬ 
tance  product  of  HyperSAS  denoted  as  Rrs^  and 
that  of  SeaPRISM  denoted  as  RteLj.  exhibits  a  spec¬ 
tral  average  AAPD  value  of  lu.2%  and  very  a 
strong  correlation,  with  an  value  of  0.972.  On  the 
other  hand,  the  spectral  average  AAPD  value  of  the 
comparison  with  the  data  processed  with  the  cur¬ 
rent  MG  algorithm  is  12.6%,  and  the  value  is 
0.958. 

Finally,  we  tested  merits  of  the  proposed  CCNY 
BRDF  algorithm  over  the  current  operational  MG  al¬ 
gorithm  by  comparing  MODIS  data  with  in  situ  data. 
The  study  clearly  shows  that  comparisons  between 
MODIS  and  in  situ  data  that  are  processed  with 
the  CCNY  algorithm  exhibit  smaller  AAPD  values 
throughout  the  spectrum,  with  a  spectral  average 
improvement  of  3.15%.  In  addition,  stronger  correla¬ 
tion  coefficients  were  also  observed  in  the  compari¬ 
son  carried  out  with  Rrs°(A)  data  processed  with 
the  CCNY  algorithm. 

In  conclusion,  we  propose  a  remote  sensing  reflec¬ 
tance  model  specifically  tuned  with  coastal  water 
lOPs  for  bidirectional  correction.  The  resulting  algo¬ 
rithm  shows  advantages  over  the  cxurent  MG  algo¬ 
rithm  on  the  tests  performed  with  simulated 
datasets,  in  situ  measured  data  sets  (for  a  one-year 
period),  and  satellite  data.  We  propose  the  use  of  this 
CCNY  BRDF  correction  algorithm  for  better  ocean 
color  data  retrievals  as  well  as  validation  and  cali¬ 
bration  for  coastal  waters. 
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